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ABSTRACT 

The radiation spectra of many of the brightest ultraluminous X-ray sources (ULXs) are dominated 
by a hard power law component, Ukely powered by a hot, optically thin corona that Comptonizes 
soft seed photons emitted from a cool, optically thick black hole accretion disk. Before its dissipation 
and subsequent conversion into coronal photon power, the randomized gravitational binding energy 
responsible for powering ULX phenomena must separate from the mass of its origin by a means other 
than, and quicker than, electron scattering-mediated radiative diffusion. Therefore, the release of 
accretion power in ULXs is not necessarily subject to Eddington-limitcd photon trapping, as long as 
it occurs in a corona. 

Motivated by these basic considerations, we present a model of ULXs powered by geometrically thin 
accretion onto stellar mass black holes. In the region closest to the hole (region I), where the majority 
of the binding energy is released, cool thermal disk radiation is Comptonized by an adjacent corona 
covering the entire surface of the disk. The amount of reprocessed thermal emission in region I is quite 
small compared to the hard coronal output since the disk behaves as a near perfect reflector of X-rays, 
a result of the intense ionizing flux which leads to an extreme state of photo-ionization. If energy 
injection takes place within an optically thin corona, the conversion of binding energy into a wind 
is hampered by Compton drag and the wind's low optical depth. Furthermore, if the magnetic field 
geometry of the corona is primarily closed, then magnetic fields of modest strength can in principle, 
prevent the launching of a wind. In the outer regions (region II), where the albedo is somewhat lower, 
thermal emission resulting from a combination of viscous dissipation within the body of the disk and 
reprocessed coronal power is emitted at relatively low temperatures, due to the large surface area. 
Within the context of the current black hole X-ray binary paradigm, our ULX model may be viewed 
as an extension of the very high state observed in Galactic sources. 

Subject headings: accretion, accretion disks - black hole physics - X-rays: binaries - X-rays: galaxies 



1. OVERVIEW OF OBSERVATIONS, CURRENT 
INTERPRETATIONS, AND PLAN OF THIS WORK 

Per detected photon, ultraluminous X-ray sources 
(ULXs) garner a disproportionate amount of attention. 
The reason for such relative popularity results from the 
two most common interpretations of their observed be- 
havior: 1) super-Eddington luminosities resulting from 
accretion onto stellar mass black holes 2) sub-Eddington 
accretion onto intermediate mass black holes. Both ex- 
planations are problematic from a theoretical point of 
view. With respect to super-Eddington luminosities, it 
is not clear how a radiatively efficient flow can release 
binding energy at a rate above the Eddington limit in 
a steady-state manner. As for intermediate mass black 
holes, a well-established evolutionary path that accounts 
for their birth docs not exist. 

Below, we summarize some important observational 
features of ULXs and previous theoretical models of their 
behavior in hopes to motivate our own work. 

1.1. Observations 

We are primarily concerned with ULXs which display 
isotropic X-ray luminosities Lx ^ 10^°ergs^^ - of order 
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10 times the Eddington limit for a \QMq black hole. The 
radiation spectra of these sources measured by Chandra 
and XMM can be roughly fit by a power law with spec- 
tral index F 2 (see Miller & Colbert 2004 for a com- 
prehensive review). For highly studied ULX sources, fits 
to the data allow for the presence of a soft thermal com- 
ponent with kT ~ 0.1 — 0.5 keV in addition to a F ~ 2 
power law. The presence of such an anomalously cool 
thermal component in many ULXs, along with the large 
inferred X-ray luminosity, provides the best argument 
for the presence of intermediate mass black holes. That 
is, large luminosities emitted thermally at relatively low 
temperatures necessarily requires a relatively large emit- 
ting area and therefore, more massive black holes. 

Table 1 lists the salient observational features of highly 
studied ULXs, selected on the basis of their large lumi- 
nosities as well as the ability of other workers to fit the 
data, in a statistically meaningfully way, with a two- 
component model consisting of a power law and soft 
multi-color disk black body. Clearly, the power law com- 
ponent in each individual source is a crucial component 
of the total X-ray power. To illustrate this point. Fig- 
ure ^ provides best-fit theoretical spectra for two sources 
listed in Table 1 . By extrapolating the power law compo- 
nent to 100 keV, the potential domination of the power 
law spectral component is further emphasized. 

1.2. Previous Interpretations 
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TABLE 1 

Spectral Properties of Selected ULXs 



Source 


Date 




r 


-Ppl/-^total 


Lx 




Reference 






(kcV) 




(10*0 erg s-i) 


(10^1 cm"2) 




NGC 1313 X-1 


2000 Oct 17 


0.23 ± 0.02 


1.76 ±0.07 


0.74=' 


0.6 ±0.1=1 


3.1 ±0.3 


1 


NGC 1313 X-2 


2000 Oct 17 


16+°-^'^ 


9 Q+0.2 
^•'-'-0.1 


0.63'= 


U.DD_Q 20 


q+3 


2 


M81 X-9 


2002 Apr 10 


26+° °2 
"■■^"-0.05 


1.73 ± 0.08 


0.85=^ 


1 1+0. 3a 
^■^-0.1 


2.3 ±0.3 


1 


M81 X-9 


2002 Apr 16 


0.21 ± 0.04 


1.86 ±0.06 


0.84=1 


1 q+0.3a 
^•'-'-0.2 


2.9 ±0.3 


1 


NGC 4038/4039 X-11 


2002 Jan 8 


0.15 ± 0.02 


1.9 ±0.2 


0.61=1 


9 1+0. 7a 
^•^-1.1 


Q+0.8 


3 


NGC 4038/4039 X-16 


2002 Jan 8 


0.19 ± 0.05 


1.4 ±0.2 


0.33=1 


^•°-i.o 


1.5 ±1.0 


3 


NGC 4038/4039 X-44 


2002 Jan 8 


15+0-02 


9 9+0.1 
^•^-0.4 


0.77=1 


1 n+^'^=i 


1 4+2-0 


3 


NGC 4559 X-7 


2003 May 27 


0.148 ± 0.006 


9 90 + 0.05 
^■^■^-0.04 


0.63<= 


2.2<^ 


5 1+1-4 


4 


NGC 4559 X-7 


2001 Jun 4 


0.12 ± 0.006 


1.8 ±0.08 


0.69<= 


1.9<= 


3.6+n 


4 


NGC 4559 X-7 


2002 Mar 14 


0.12 ± 0.01 


2.13 ±0.08 


0.54<= 


3.2<^ 


5.7t°i 


4 


M74 X-1 


2001 Oct 19 


0.31 ± 0.05 


l-o3_Q 5g 


0.78=1 


0.4=1 


0.48 


5 


Holmberg II X-1 


2002 Apr 10 


141+0.018 




0.83^ 


2.0=1 


1 6+°-2 


6 


Holmberg II X-1 


2002 Apr 16 


n 1 90+0.022 




0.74d 


1.7=1 


1.4 ±0.3 


6 


Holmberg II X-1 


2002 Sep 18 


120+°°^^ 


^■"='-0.08 


0.68^ 


0.5=1 


1 4+0-5 


6 


Holmberg II X-1 


2004 Apr 15 


0.20 ±-1-0.02 


2.64 ±0.03 


0.91^ 


1.2'* 




7 


NGC 5204 X-1 


2003 Jan 6 


0.21 ±-1-0.03 


1.97 ±0.07 


0.82=i'f 


0.44=1 


0.781°;- 


8 


NGC 5204 X-1 


2003 Apr 25 


^^•^^-0.03 




0.69=i'f 


0.55=1 


1 99+0.34 

-l-^^-0.29 


8 


MIDI XMM-1 (P13/H19) 


2002 Jun 4 




1 41+0.25 
^■^^-0.18 


0.74=1 


0.27=1 




9 


NGC 7771 X-2 


2002 June 21 


16+°-l2 
^•^"-0.03 


^■°'-0.26 


0.80=1 


3-86±^;2i" 


3.18l^;« 


10 



References. — (1) Miller et al. (2004a); (2) Miller ct al. (2003); (3) Miller et al. (2004b) ; (4) Cropper ot al. 2004; (5) Krauss ct al. 
2005; (6) Dewangan et al. 2004; (7) Goad ct al. 2005; (8) Roberts al. 2005; (9) Jenkins ct al. 2004; (10) Jenkins et al. 2005 



=iCalculatcd over the energy range 0.3-10 keV. 

''Calculated over the energy range 0.2-10 keV. 

•^Calculated over the energy range 0.3-12 keV. 

"iCalculated over the energy range 0.3-2 keV. 

° Their definition of F is equal to our definition minus one 

^We calculated the ratio based on parameters of their best-fit model. 



Several models have been put forth to explain ULX 
phenomenon. For Lx }^ lO"'" erg s^^, accretion onto in- 
termediate mass black holes (IMBHs) receives much at- 
tention. If interpreted as emission from an optically thick 
accretion disk extending to the inner-most stable circular 
orbit of the black hole, the temperature of the soft ther- 
mal component inferred from the spectra of these sources 
are consistent with black hole masses of ^ 10^ — IO'^Mq. 
However, the large fraction of non-thermal emission in 
these sources (see Table 1) brings into question the ro- 
bustness of this interpretation. That is, the IMBH ar- 
gument is straightforward only when the cool thermal 
component is associated with the majority of the bolo- 
metric luminosity, rather than a sub-dominant fraction 
of the gravitational power. 

Alternatively, ULXs may result from accretion onto 
~ IOM0 black holes. In this case, the emission must 
either be beamed so that the intrinsic energy release 
is less than that inferred (King 2001) or the Edding- 
ton limit must be exceeded. Begelman (2002) suggests 
that the non- linear development of the photon bubble in- 
stability (Gammie 1998; Blaes & Socrates 2003; Turner 
ct al. 2005) might produce channels in the accretion 
flow, allowing the disk luminosity to exceed the Edding- 
ton limit since radiative trapping is overcome. In this 
work we propose an alternative scenario for generating 
steady state super-Eddington luminosities from stellar 
mass black holes, with the intent of reproducing the 



hard, non-thermal emission that dominates the spectra 
of many bright ULXs. 

1.3. Plan of this work 

In the light of the rough observational guidelines out- 
lined in iil.ll for ULXs, an accretion model describing 
ULX behavior must contain certain indispensable in- 
gredients. Most importantly, the observed large fluxes 
should be released non-explosively. A good candidate 
is radiatively efficient accretion onto a compact object, 
whose birth and origin are well understood. Such an 
accretion flow must account for the large apparent coro- 
nal output while simultaneously accounting for the sub- 
dominant and relatively cool 0.1-0.5 keV thermal emis- 
sion component. 

In the next section we provide the framework in which 
the above observational constrains are satisfied. In 30 
the central idea of this work is motivated and argued. 
That is, the release of accretion power resulting from a 
two-phased accretion flow is not bound by the Edding- 
ton limit since radiative diffusion, mediated by electron 
scattering, is not primarily responsible for the removal 
of randomized binding energy from the mass of its ori- 
gin. The feasibility of confining coronal magnetic fields 
as well as the effects of Compton drag on the outgoing 
optically thin wind are discussed in ^2.2\ and H2.3I re- 
spectively. The formation of ULX spectra in terms of a 
corona+disk model is considered in ^and we conclude 
in g| 
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2. THE BASIC IDEA AND INITIAL ESTIMATES 

Perhaps the most important assumption of one-zone 
models of accretion is that the randomized gravitational 
binding energy is dissipated locally with the mass of its 
origin. In the case of thin accretion disks, the gravita- 
tional power is transported vertically via radiative diffu- 
sion from the midplane of disk, where most of the mass 
resides. As a result, super-Eddington accretion rates do 
not proportionally lead to super-Eddington luminosities 
since the radial inflow time becomes short relative to the 
time it takes for a photon to escape from the disk mid- 
plane. That is, the photon trapping radius moves out- 
ward with accretion rate, negating the increase in avail- 
able binding energy with a decrease in the overall radia- 
tive efficiency. 

Another consequence of the local dissipation assump- 
tion is that the emergent radiation spectrum will be ther- 
mal and approach the black body limit. However, spec- 
tral components that probe the central engines of ac- 
tive galactic nuclei (AGN) and X-ray binaries (XRBs) 
indicate the presence of powerful non-thermal coronal 
emission accompanied by a soft quasi-thermal compo- 
nent. A common interpretation of AGN and XRB spec- 
tra is two-phased accretion onto a black hole (Haardt & 
Maraschi 1991, 1993, hereafter HM91, HM93; Svennsson 
& Zdziarski 1994, hereafter SZ94). In these models, a 
fraction / of the gravitational energy release is assumed 
to dissipate in a hot diffuse corona, above the main body 
of the disk, away from the majority of the accreted mass. 
It is typically thought that the viscous stress, assumed to 
be magnetic in nature (Balbus & Hawley 1991, hereafter 
BH91), transports angular momentum and initially ran- 
domizes the gravitational binding energy near the mid- 
plane. The magnetized fluid elements, which are buoyant 




Energy (keV) 

Fig. 1. — The best fit MCD + power-law model spectra are 
plotted for NGC 4038/4039 X-44 (blue) and M89 X-9 (2002 Apr 10; 
red). The curves represent the flux in the unabsorbed total model 
(solid), MCD component (dashed), and power-law component (dot- 
dashed) for the parameters listed in Table 1. Note that the models 
were fit with a power-law component over the 0.3-10 kcV bands, 
but we have extrapolated the models assuming 100 keV cut-off to 
the power-law for the purpose of illustration. We have also rescaled 
the overall flux of M89 X-9 for illustrative purposes as well. 



with respect to their surroundings, dissipate - perhaps 
through magnetic reconnection or shocks - above the disk 
where the reconnection time and Mach number are large. 

For large values of /, often necessary to fit XRB and 
AGN data, the non-local dissipation of binding energy 
can have profound consequences on the disk structure 
as well as the net radiative efficiency. By introducing a 
mechanism of energy transfer other than radiative diffu- 
sion, the relevance of the Eddington limit is diminished 
(but not eliminated entirely). Below, we qualitatively 
discuss why super-Eddington luminosities may result for 
two-phased accretion onto a black hole. 

2.1. Altered disk structure and trapping radius 

Following SZ94, we assume that all of the accretion 
takes place in the main body of the disk. Since we are 
interested in explaining ULX behavior with stellar mass 
black holes, we restrict our discussion to accretion rates 
above the Eddington limit for a 3-20 M0 black hole, 
implying a radiation pressure dominated disk. 

The difference in structure between a standard thin 
disk and that of a disk -|- corona is due the condition of 
radiative equilibrium 

^. = .T,V-(i-/)Q-(i-/)|^^| (1) 

where Q is the viscous dissipation rate per unit area. All 
changes to disk quantities due to coronal dissipation de- 
pends on (1 — /) and for reference, their values (as well 
as definitions of the various disk quantities) with respect 
to standard thin disk scalings are given in Appendix 
1X1 Most notably, when compared to its classic one-zone 
counterpart, the disk scale height Hd decreases by a fac- 
tor (1 — /) and therefore, the cool disk remains thin as 
long as an increase in m is proportionately matched by an 
increase in /. Furthermore, decreasing the internal dis- 
sipation increases the disk midplane density increases by 
a factor (1 — f)~^. This has the important consequence 
that now, the radial inflow velocity becomes 

VB= , ^„ = voB.il -ff'^M'il-ff (2) 
47r RHdPd 

where wofl is the classic thin disk value (SS73). As long 
as increasing m is matched by an increase in /, then 
remains fixed, implying a reduction in the level of radia- 
tive trapping. To quantify this statement even further, 
we define the trapping radius Rtr, defined as the radius 
at which the radial inflow time is equal to the vertical 
diffusion time, 

Rtr^^^RAl-ff'^m. (3) 

As expected, an increase in m is matched by an increase 
in Rtr, directly leading to a decrease in radiative effi- 
ciency. For / ~ at the Eddington accretion rate, the 
trapping radius Rtr ^ lORg for a Schwarzschild hole. 
This implies that thin disk accretion cannot maintain 
its radiative efficiency beyond the Eddington accretion 
rate if the release of binding energy is governed by elec- 
tron scattering-mediated radiative diffusion. However, 
for large values of / the binding energy of the gas is 
removed from large to small optical depths by a mech- 
anism other than radiative diffusion, perhaps magnetic 
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buoyancy or waves. In this case, the modified radial in- 
flow time R/vR must be compared to the characteristic 
timescale of say, buoyant transport. 

Without a coherent theory of relativistic accretion disk 
turbulence, it is difficult to estimate what the character- 
istic speed of magnetized buoyant transport vb will be. 
A good guess for the upper limit oivs is the sound speed 
of the disk at the midplane Cg. Let us define an inward 
advection time Tadv = R/vr and a timescale for upward 
buoyant transport tb = Hd/cs. As long as the ratio 
Tadv/^B > 1, the accretion power of the fiow stands a 
chance of escaping before being consumed by the black 
hole. In terms of disk quantities, we have 



TB ^ rh am?{l-ff 

If we (somewhat arbitrarily) take a to be ^ 0.1, then the 
above ratio assumes a minimum value of Tadv/TB ~ 5 at 
r ~ 12, the radius of maximum light, for an accretion rate 
rate of m = 175 - ten times the Eddington rate - for a 
choice of (1 — /) ~ 0.1. Therefore, the efficiency of the 
disk remains unchanged as long as an increase in the ac- 
cretion rate rn is proportionately matched by a decrease 
in (1 — /). Even if the rate at which randomized gravita- 
tional binding energy leaks upwards decreases by a factor 
of 5, the disk may maintain its radiative efficiency for 
accretion rates of order 10 times the Eddington rate. 

2.2. Magnetic field requirements 

Any viable accretion model that steadily generates 
super-Eddington luminosities must fulfill two require- 
ments. First, the accretion flow must remove its ran- 
domized binding energy before being advected into the 
hole in order to maintain a large radiative efficiency. Fur- 
thermore, the majority of that binding energy must be 
converted into radiation rather than a mechanical out- 
flow. The central point of this paper, encapsulated by eq. 
(@J, addresses the first requirement. That is, magnetic 
buoyancy is an ideal mechanism for the non-dissipative 
removal of binding energy from the bulk of the flow. Ad- 
dressing the second requirement is conceptually more dif- 
ficult and again, we appeal to the existence of magnetic 
fields. 

Irrespective of the geometry of the emitting surface, 
a super-Eddington radiation flux produces an outward 
opposing acceleration that is larger than the attrac- 
tive inward acceleration provided by gravity (see e.g., 
Abramowicz et al. 1980). Therefore, an additional force 
or stress is required in order for our super-Eddington ac- 
cretion model to maintain its radiative efficiency. Begel- 
man (2002) suggests that strong magnetic fields might 
prevent surface regions of low optical depth from being 
blown off by a super-Eddington fiux generated near the 
disk midplane. Furthermore, several authors have noted 
that accretion flows which liberate a sizable fraction of 
their binding energy in a Comptonzing corona inevitably 
require modest magnetic fields in order to mitigate the 
escape of pairs by the act of confinement (Svensson 1984; 
Zdziarski 1985; White & Lightman 1989). We examine 
a similar scenario in the context of the super-Eddington 
flux arising from energy injection in an optically thin 
corona. If indeed the conflnement is provided for by 
magnetic fields, then there are some straightforward con- 
straints regarding its topology and strength. 



Charged particles are free to move along magnetic field 
lines. As long as the held lines are anchored in the 
disk, the supposed site of their generation, their magnetic 
tension is capable of confining fluid elements. Numeri- 
cal simulations of stratified accretion flows indicate that 
the coronal field geomtry is primarily closed and toroidal 
(Miller & Stone 2000; Hirose et al. 2004). However, it 
is not clear whether or not field lines are anchored deep 
in the disk. To gain an order of magnitude estimate of 
the acceleration a confining magnetic field can impart to 
a electrically conducting coronal fluid element, we define 
gM to be a characteristic magnetic acceleration as 




which should be compared to the magnitude of the grav- 
itational acceleration 

and the radiative acceleration 

.... = ^F.^2xlO-(-)(-).-3_ 

which has a maximum value of grad 2 x 10^'^ for 
m = 175atr~12. Note, m is the accretion rate scaled 
to the Eddington rate Msdd = ^irGMmp/ {cut) and m is 
the black hole mass in units of Mq . If the radiative force 
in the corona is balanced by gradients in the magnetic 
field, then the ratio gM 1 9rad > 1 and the quantity grad/ 9 
is approximately the factor by which the Eddington limit 
is surpassed. If we assume that He and R scale in pro- 
portion to the mass of the hole, then like the charac- 
teristic magnetic and radiative acceleration gM and grad 
are oc Afgj^ and therefore in principle, similar Eddington 
ratios can be achieved for supermassive black holes. 

The magnetic field strength or pressure in the corona 
can be constrained as well. First, the magnetic pressure 
of the corona Pm.c ~ B1/%tt must be at least compa- 
rable to the coronal pressure. The other requirement is 
that the coronal magnetic pressure can not exceed the 
magnetic pressure at the disk midplane. Following the 
scalings of SZ94, the required coronal field strength is 
given by equating Pra,c to Prad,c (see Appcndi jX)) . which 
yields 

1 /2 

> 2.4 X 10^ mi/2Ai/2^-3/2 q (g) 

for / 1. By assuming a magnetic origin for the accre- 
tion torque, an estimate of the magnetic field strength in 
the disk midplane is obtained by setting Pm,d ~ trjP ^ 
aPd, giving us 

Bd - 1.6 X 10^ (1 - /)-i/2 ^'^ G. (9) 

For m = 175, / = 0.9, and A ^ 1 at a radius of r ~ 12 
the ratio Bd/B^ ~ 3, which is its lowest value. Thus 
it seems possible that if only a fraction of the dynamo- 
generated flux resides in the corona, the plasma there 
can be contained despite the upward force exerted upon 
it by the super-Eddington radiation field. 

As a check of internal consistency, the coronal magnetic 
field given by eq. lO must lead to an Alfvcn speed va 
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large enough to satisfy cq. jSJ. In order to evaluate 
VA, knowledge of the coronal density is required. For 
a radiation pressure dominated corona, a characteristic 
density is roughly given in Appendix ^ leading to an 
Alfven speed of ~ c for the choice of disk and coronal 
parameters previously mentioned. 

2.2.1. Hierarchy of stresses: Or, why coronae are 
essential for super- Eddington luminosities 

Consider the following relevant thought experiment. 
Imagine injecting thermal energy into the Sun at a rate 
of ~ 10"^^ ergs~^, which is ~ 10 times its Eddington rate. 
Regardless of the depth at which this enormous deposi- 
tion of energy takes place, the local pressure will simi- 
larly increase by a factor of ^ 10. This can be under- 
stood by noting that the ratio of radiation to gas pressure 
Prad/Pg is roughly a constant throughout the star and 
approximately equal to its Eddington factor such that 

Prad/Pg ~ L^/LEdd- 

If the injection takes place uniformly at a great depth 
near the core, then the resulting radiation stress is ^ 10 
times larger than any other stress in the entire star. 
Thus, the bulk of the star has no choice other than 
to move outward and form a massive wind that carries 
away the injected energy. However, if the energy injec- 
tion uniformly takes place at the surface and near the 
photosphere, then the situation can potentially, be quite 
different. That is, a radiation pressure supplying a flux 
~ 10 times the Eddington limit at the surface of the 
Sun is only ~ 10^^^ that of the gas pressure in the core. 
Hence, super-Eddington energy injection at the surface 
does not appreciably change the structure of the star. If 
the stresses within the star re-arrange themselves only 
slightly, then a super-Eddington flux due to energy in- 
jection at the surface can easily be confined as to sti- 
fle the production of a non-radiative super-Eddington 
wind. To further clarify this point, take the example of a 
sunspot, whose magnetic stresses are supposedly gener- 
ated by highly subsonic (and highly sub-virial) turbulent 
and shearing motions at the base of the solar convection 
zone. For typical sunspot field strengths ^ a few xlO'' 
G, closed field lines at the photosphere could in principle, 
confine matter that is accelerated by radiation pressure 
of order ~ 10 times Eddington. Despite the fact that 
only an infinitesimal amount of the Sun's internal energy 
resides in its convection zone, subsonic motions in these 
regions may lead to stresses large enough to contain re- 
gions of the photosphere experiencing a super-Eddington 
radiation pressure. 

The situation described in the previous section with 
respect to ULXs powered by super-Eddington accretion 
disks is quite similar to the solar thought experiment 
outlined above. Even for super-Eddington luminosities, 
the ratio of coronal to midplane pressure PcjPd is quite 
small. In the disk case, the magnetic accretion stress 
aPd originates from motions that are subsonic in the 
midplane by a different mechanism than, but analogous 
to, the generation of field loops in the shear layer at 
the base of the solar convection zone. The arguments 
surrounding eqs. and jnj indicate that if even a small 
fraction of the magnetic energy leaks into photosphere, 
magnetic confinement is a plausible mechanism for the 
suppression of a super-Eddington wind. 



2.3. Compton- driven wind 

In the absence of closed magnetic field lines, it is com- 
monly thought that the majority of the gravitational 
power for a super-Eddington accretion flow either is ad- 
verted into the black hole or is converted into mechan- 
ical energy that drives an outflow. In this two-phased 
accretion model, it is unlikely that a significant amount 
of the binding energy is advcctcd into the hole for suffi- 
ciently large values of /, the fraction of power released 
in the corona. Therefore, most of the radiative luminos- 
ity would be expected to be converted into the kinetic 
power of the outflow. However, this conclusion follows 
from the assumption that the majority of the energy in- 
jection takes place beneath the photosphere. 

The properties of continuum driven winds follow from 
conservation of momentum (Lamers & Cassineli 1999) 

and mass 

M = A.-kR^pv. (11) 

Here, is the number density of the ions and n+ is 
the number density of pairs. For the sake of simplic- 
ity, we have assumed spherical symmetry. Since we are 
considering super-Eddington luminosities, the forces due 
to gravity and gas pressure can be neglected outside the 
sonic point. Integrating over space, eq. H10|l then be- 
comes 

Mwo, ~ Tp- f 1 + 2^ ) , (12) 

where Tp and Woo are the plasma optical depth and out- 
flow velocity at infinity, respectively. Note that the ra- 
tio n+jrip is approximated as constant. The above ex- 
pressions are accurate for non-relativistic velocities. Due 
to the large radiative acceleration, a relativistic outflow 
seems inevitable. However, Compton drag limits the 
Lorentz factor to moderate values, such that F^ > 1 
(Madau & Thompson 2000). 

For Vao ~ c, the kinetic luminosity in the wind is 
simply 

-TpL (^1 + 2^^ . (13) 

Thus, the kinetic power resulting from super-Eddington 
energy injection is directly proportional to the total 
plasma -I- pair scattering optical depth and in principle, 
may be significantly smaller than the photon luminosity. 
These arguments are simple. In order to truly quantify 
the rate of loss of mechanical energy, the spatial distri- 
bution of energy injection must be specified, potentially 
allowing for a detailed calculation of the coronal struc- 
ture, pair equilibrium, and resulting spectrum. 

3. CORONAL STRUCTURE AND ULX SPECTRA 

The radiation spectra of high luminosity ULXs are of- 
ten dominated by a relatively flat T r-j 2 power law. Also, 
the fractional thermal emission, responsible for ~ ^ to 
of the luminosity, is typically modeled by a black body 
or multi-temperature disk black body with a character- 
istic temperature of 0.1-0.5 keV. Along with producing 




REGION I ' ^ REGION II 



Fig. 2. — Sketch of our stellar mass black hole ULX accretion model. In both Regions I and II, gravitational binding energy is removed 
from the disk in the vertical direction before it is converted into heat. However, accretion power may dissipate locally in Region II without 
being advected into the hole. The majority of the gravitational power is released in Region I primarily in the form of hard Comptonized 
photons. Throughout the flow, cool photons from the optically thick disk are Comptonized in the surrounding hot diffuse corona. Hard 
X-ray photons incident upon the disk are nearly perfectly reflected in Region I while in Region II, they are photo-electrically absorbed 
by the disk. In the corona, the production of an energetically dominant super-Eddington wind is hindered by the presence of confining 
magnetic fields, a low optical depth, and the effects of Compton drag. Note that this figure is not to scale. 



a luminosity of ^ 10^° — lO^^ergs"^. these two obser- 
vational constraints must be reproduced by an accretion 
flow spectral model that describes ULX behavior. 

Our model consists of a two-phased accretion flow di- 
vided into two distinct geometrically thin regions (see 
Figure [^. The inner portion of the disk, or region I. 
is responsible for liberating the majority of the system's 
gravitational binding energy, while the outer portion, re- 
gion II, powers the sub-dominant thermal spectral com- 
ponent. 

In region I, we adopt the slab geometry for the 
disk-|-corona i.e., the simplest configuration possible. In 
order to reproduce the spectra of luminous ULXs, re- 
gion I cannot be a prodigious site of soft thermal photon 
power. Therefore, the emergent spectra from the region 
of deepest gravitational potential must resemble a F 2 
Comptonized power law with only a small fraction of the 
luminosity in the form of thermal emission such that ra- 
tio of hard coronal to soft power Lc/Ls ~ 10 in region 
I. 

As described in the previous section, if randomized 



binding energy is deposited in a corona in excess of the 
Eddington rate, hydrostatic balance cannot be achieved 
unless forces that are non-gravitational in nature, such 
as magnetic stresses, are present. In the absence of mag- 
netic fields, the corona becomes geometrically thick and 
a wind develops. Calculating the structure of the corona 
by accurately including the relevant processes that ac- 
company magnetic field generation and transport, non- 
equilibrium radiation physics, and winds is an impossibly 
difficult problem, with no clear way as to how to pro- 
ceed. However, some simplification is achieved by real- 
izing that Comptonization of ambient soft photons must 
be the radiative mechanism which mediates the release 
of energy in the corona. That is, relatively few physical 
parameters, such as optical depth and temperature, are 
required to reproduce the observed spectral components 
for a given coronal geometry. Though such an approach 
lacks almost any sense of predictive power, it provides a 
rough outline as to what the flow geometry and energet- 
ics might be. 

In what follows, we briefly discuss the composite phys- 
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ical requirements of the accretion flow necessary to re- 
produce the observed spectra. 

3.1. Spectral index, coronal temperature, and optical 
depth of the inner region 

An almost universal feature of the spectral energy dis- 
tributions (SEDs) of the relativistic accretion flows that 
power active galactic nuclei (AGN) and XRBs is a pre- 
dominantly flat broadband X-ray power law which results 
from thermal Comptonization of hot coronal electrons. 
In AGN and black hole XRBs in the low/hard state, the 
spectral index F < 2, with greater variation seen in black 
hole XRBs. 

In the case of low/hard state spectra of galactic black 
hole candidates such as Cyg X-1, both the local and 
disk- integrated values of Lc/Lg ^ 1. The most pop- 
ular disk-|-corona geometry that is invoked to explain 
these sources is a truncated disk -I- hot inner spherical 
corona model (Poutanen, Krolik, & Ryde 1997; Esin et 
al. 1998). The inner hot regions are thought to resem- 
ble the ADAF/CDAF accretion models that are thought 
spontaneously occur at low accretion rates (Narayan & 
Yi 1995; Quataert & Gruzinov 1999). However, there 
are other explanations for the low/hard state that can- 
not be ruled out and both involve altering the patchy 
corona model of IIMG93 by reducing the soft reprocessed 
disk emission. If the overwhelming majority of the flow's 
gravitational power is liberated in relatively small active 
coronal regions, then the ionizing flux incident upon the 
disk will be so large that the atomic X-ray absorptivity of 
the disk becomes vanishingiy small and disk acts almost 
as a perfect reflector (Ross & Fabian 1993; Ross, Fabian. 
& Young 1999; Ballantyne, Ross, & Fabian 2001). An- 
other mechanism that reduces reprocessed disk emission 
in the patchy corona model is outward vertical motion 
of the active regions at sub-relativistic velocities (Be- 
loborodov 1999). 

In region I of our ULX model, the local and disk- 
integrated value of Lc/Lg 1 as in the case of black 
hole XRBs in the low/hard state. For a spectral index 
r ~ 2, equal photon power is emitted at all photon en- 
ergies bounded by the input seed photon energy and the 
coronal electron temperature. Therefore, the majority 
of the disk photons interact with the corona only a few 
times, if any, and the relatively few photons that stay in 
the corona the longest carry away with them most of the 
radiated power. The likelihood for a photon to escape 
the corona, which varies inversely with the optical depth 
T, is appropriately balanced by the average photon en- 
ergy shift per scattering AE/E ~ AkBTc/rne(? in order 
to reproduce a given spectral index F (see e.g., Fermi 
1949 within the context of cosmic ray acceleration in the 
interstellar medium). Thus, if the ksTe ~ 100 keV, then 
only of the input photons would have to leave the 
system at that energy for a soft seed photon energy of 
~ 1 keV and an outgoing spectral index F ^ 2. Such 
harmony can be achieved for r ^ 1 since /S.E/E < 1 if 
ksTe ~ 100 keV. That is, a F ~ 2 Comptonized spec- 
trum is viable only when the propensity of escape per 
scatter is proportionally countered by a gain in energy 
for every scattering event. 

For flat F ^ 2 Comptonized spectra, the amplification 
Lc/Ls can be estimated by calculating the number of 
scattering orders Ns required to fill out the Comptonized 



power law, since the location of each scattering order 
in the SED carries equal amounts of energy as the soft 
input. By choosing soft seed photon energies of order ~ 1 
keV, roughly an Eddington's worth of power is emitted 
from the cold optically thick disk. From the argument 
above, F ~ 2 for coronal temperatures of order 100 
keV and an optical depth 1. For these parameters, we 
obtain iV^ ~ Lc/Lg < 10. Thus, an input soft seed X- 
ray source typical for optically thick Eddington-limited 
accretion onto IOMq black hole is amplified up to ^ 10 
times its intrinsic thermal power if there is an adjoining 
corona with ^ 100 keV and t ^ 1 covering the entire 
disk. 

The somewhat prosaic arguments given above are in 
rough agreement with detailed Comptonization calcula- 
tions of other authors (Pietrini and Krolik 1995; Stern et 
al. 1996). For example, Pietrini and Krolik (1995) de- 
duced a scaling relationship between T^, r, Lc/La, and 
F that reproduce the desired spectral properties of our 
ULX coronal model. Similar to the low/hard state of 
black hole XRBs, our ULX model requires a low level of 
soft photon input. We cannot appeal to a truncated disk 
-f ADAF/CDAF scenario because large accretion rates 
and high radiative efficiencies are necessary for ULXs. 
In order to reduce the amount of reprocessed thermal 
power, the albedo of the disk must approach unity due 
to the intense super-Eddington ionizing fiux. In what 
follows, we discuss the validity of this assertion. 

3.2. Albedo 

In order to reproduce spectra dominated by a Comp- 
tonized power law in the slab geometry, the Comptonized 
photons incident upon the disk must avoid photo-electric 
absorption. If a significant fraction of the hard X-ray flux 
is absorbed by the disk, the reprocessed thermal emis- 
sion can in principle be comparable to the local X-ray 
flux such that L^jLs ^ 1 in the slab geometry (Haardt 
& Maraschi 1991). Thus, either the disk albedo must 
approach unity as a result of photo-ionization (Ross & 
Fabian 1993) or the coronal electrons flow away from 
the disk at sub-relativistic to relativistic velocities (Be- 
loborodov 1999). 

In a super-Eddington model of ULXs, the ionizing 
coronal flux is intense in comparison to say, ionizing disk 
models of black hole XRBs in the low/hard state. An 
increase in disk albedo implies a change in the disk's role 
from that of an absorber and re-emitter to a reflecting 
X-ray mirror. If the disk albedo Od approaches its limit- 
ing value such that = 1, nearly all of the metal ions 
in surface of the disk must be fully ionized. 

In Appendix ^ Od is shown to closely approach unity 
near the inner-edge of the disk for m =175 while sharply 
dropping off at larger radii. By r ~ 60, outside of which 
1/5 of the gravitational power is released, Od ~ 1/2. 
Therefore, one expects a signiflcant amount, at least 
compared to the Eddington limit, of thermal emission 
to be emitted from ULX flows, despite the majority of 
the flow's gravitational power being released via Comp- 
tonization. 

Eq. (|B5p raises an interesting point. For accretion 
disks that power AGN, the abundant CNO metals are 
not fully ionized since the thermal disk temperatures are 
significantly lower. As a result, Az ~ 10^^ — 10^'^, imply- 
ing that the disk cannot act as a perfect reflector and the 
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spectra will be different than the XRB case in that the 
reprocessed thermal spectral component is more likely to 
be comparable in magnitude to the Comptonized power 
law. 

3.3. Thermal component 

The inferred thermal component of bright ULXs are 
often cited as the most compelling evidence for the exis- 
tence of intermediate mass black holes due to their rel- 
atively low temperatures. Roughly speaking, the most 
luminous sources can be fit just as well with a simple 
power law rather than a sub-dominant disk -f dominant 
power law component. Therefore, explaining the phys- 
ical processes responsible for the inferred soft thermal 
component is not central to our understanding of ULXs. 
Nevertheless, we describe the manner in which soft ther- 
mal emission may emerge from our ULX model. 

In our model, thermal disk emission primarily origi- 
nates from region II. Take the example of a lOM© black 
hole accretion flow which resembles our ULX model ac- 
creting at 10 times the Eddington rate rriEdd- For a non- 
rotating hole, the brightest region of the flow is located 
at r ^ 12 and since gravitational power falls of as oc r~^, 
roughly 90% of the accretion power is liberated within 
r ~ 120. Assume that this luminous inner portion of the 
flow constitutes region I. Outside this radius in region II, 
further assume that the dissipation and release of grav- 
itational binding energy proceeds in the classical SS73 
sense such that an Eddington's worth of gravitational 
power is emitted there as a result of viscous dissipation. 
In this case, the disk remains thin at r ^ 120 such that 
H / R < 1 , the diffusion time is short compared to the in- 
flow time, and the radiation spectra is roughly thermal. 
Compared to a disk with m = niEdd around the same 
black hole, the temperature of the thermal emission from 
region II in the ULX case is cooler due to its larger emit- 
ting surface. Typical inner-disk effective temperatures 
for an Eddington accretor around a IOMq black hole are 
ksTeff ~ 0.7 keV, whereas the m — lOniEdd ULX model 
radiates an equal amount of power at a lower tempera- 
ture of - 0.2 - 0.3 keV. 

An energetically dominant corona is required in order 
to generate more than an Eddington's worth of thermal 
power from region II in order to avoid the trapping and 
thickness problem. Therefore, super-Eddington values of 
thermal power must result from the reprocessing of coro- 
nal photons inside the trapping radius. If wc increase 
the accretion rate to bOmEdd, then for the same emit- 
ting area for region II discussed above, ~ 5 times an 
Eddington's worth of thermal power will be emitted at a 
temperature of ~ 0.4 — 0.5 keV outside of r 120. Evi- 
dently, our ULX model possesses the ability to produce 
over an Eddington's worth of thermal power for a IOMq 
black hole at relatively low temperatures in what we have 
termed region II. 

To mitigate the radiative inefficiencies due to the 
production of a powerful wind, which inevitably ac- 
companies super-Eddington thermal fluxes, region IPs 
corona must possess and ordered, predominantly az- 
imuthal magnetic field. By use of eqs. © and we 
see that the required coronal field strength Be in region 
II is small compared to value of the field in the midplane 
such that Be/ Bd < -h for sufficiently large values of /. 



4. DISCUSSION AND SUMMARY 

Many bright ULXs exhibit power-law emission with 
hard F < 2 spectral indexes. In fact, recent Chandra sur- 
vey data indicates that ULXs on average have F ~ 1.8 
with significant scatter, but independent of luminosity 
(Swartz et al. 2004). If ULXs are simply scaled up 
versions of Galactic black hole XRBs, then their F < 2 
power-laws are consistent with the low/hard state spec- 
tra of Galactic sources (Remillard & McClintock 2003). 
For this scenario, the accretion fiow is expected to be 
sub-Eddington such that Lx ^ 10~^ Ledd (e.g. Mac- 
carone 2003), implying black hole masses > 10 Mq. In 
the ADAF/CDAF model of the hard state, the thermal 
emission is produced by a truncated disk at relatively 
large radii and correspondingly low temperatures. If we 
adopt this picture for ULXs, then the temperature of 
the truncated disk component would be a factor of 3-7 
times lower than the ^0.1 keV truncated disk temper- 
ature found in Galactic black hole XRBs such as Cyg 
X-1 (Gierlinski et al. 1997; Poutanen, Krolik, & Ryde 
1997). As a result, the thermal components near ~ 0.1 
keV which are inferred from fits to some bright ULXs 
would not be produced in the truncated disk model. 

Thus, an explanation of the soft thermal component in- 
ferred from bright ULX spectra requires a disk -I- corona 
geometry that is different from the above picture for 
Galactic XRBs, whether or not the central object is 
an IMBH. If ULXs are indeed powered by an IMBH 
ADAF/CDAF- type accretion flow, then the soft spec- 
tral component must come from another source such as 
diffuse nebular emission. In any case, for hard F < 2 
power-laws, the soft component of ULXs cannot be uti- 
lized to deduce the mass of the black hole in a straight- 
forward manner, unless the overall accretion geometry is 
markedly different from that of Galactic XRBs. 

Our accretion model is primarily motivated by the fact 
that many bright ULXs are dominated by hard non- 
thermal photon power. This strongly suggests that ra- 
diative diffusion is not responsible for transporting the 
majority of the gravitational binding energy out of the 
accretion flow. Photon trapping, a mechanism which 
caps the luminosity of standard thin/shm disks, need 
not limit the flow's output to the Eddington value. For 
super-Eddington energy injection rates, the outward ra- 
diative acceleration outstrips the local gravitational ac- 
celeration. If reasonably strong magnetic stresses are 
present in the corona, the generation of a continuum- 
driven wind is prevented as long as the coronal field line 
are closed and anchored in the disk. In the absence of 
confining magnetic fields, the resulting wind inefficiently 
converts photon energy into outfiowing mechanical power 
due to the effects of Compton drag and a low optical 
depth. Therefore, the majority of the super-Eddington 
accretion power is not necessarily converted into mechan- 
ical energy and thus, the flow's radiative efficiency is left 
relatively untouched. In the corona, randomized gravi- 
tational binding energy is converted into radiative power 
via Compton amplification of soft seed photons. As a 
result of the large downward ionizing flux deep in the 
gravitational potential, the disk's albedo is sufficiently 
high such that it behaves as a reflecting X-ray mirror 
and relatively little thermal emission emerges from the 
luminous inner regions. Further away from the black 
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hole a combination of reprocessed emission and viscous 
dissipation may produce roughly an Eddington's worth 
of photon power at low thermal energies since the emit- 
ting area is comparatively large. 

Phenomenologically, our ULX model may be viewed 
as a natural extension of the commonly accepted picture 
of how Galactic black hole XRBs evolve with luminos- 
ity (e.g. Done & Gierlihski 2004). For soft state XRBs, 
the ratio of coronal to soft power Lc/Lg increases with 
luminosity. In the "high state" Lc/Lg ~ 0.1 whereas 
in the "very high state," typical values of Lc/Lg ~ 1. 
Furthermore, L ~ Q.lLEdd in the high state while the 
luminosity of the very state approaches LEdd and inter- 
estingly, Lc/Lg ^ 10 for ULXs with L ~ 10 LEdd- 

Of course, there are many theoretical uncertainties 



uniquely haunting the ULX model presented in this work. 
In order to overcome the Eddington limit by a factor 
of ~ 10, the fraction of binding energy released in the 
corona / ~ 0.9 so that the bulk of the flow can expel 
its binding energy before being advected into the hole, 
while remaining geometrically thin. If magnetic fields 
are utilized as a confining mechanism in the corona, their 
strength, geometry, and stability must be considered in 
further detail. The net radiative efficiency of the flow 
hinges upon the assertion that the majority of the ac- 
cretion power does not escape in a mechanical form. In 
order to verify this claim, the vertical structure of the 
corona resulting from the intense Comptonizing radia- 
tion field, must be be determined. 



APPENDIX 
A: DISK SCALINGS 

In this appendix, we provide disk scalings relevant for the disk midplane following the work of SZ94. When the 
accretion rate is large we choose to ignore radial advection, an assertion that is justified as long as the fraction of the 
energy dissipated in the corona / — 1 i.e., the regime relevant for ULX behavior. 

In the radiation pressure dominated limit, the midplane disk quantities are given by 

|^ = ^-^-^(l-/)^, (Al) 



8 TOr, 



-im-V3/2(l-/)- 



DM' 



Td = -a 



(1-/)" 



, CB 



(A2) 
(A3) 



Prad,d = Pd = l^^^-'r-^'i^ fr'j- (A4) 

The factors A — V enforce the no-torque inner boundary condition and take into account the effects of general relativity 
(Riffert and Herold 1995; Hubcny and Hubeny 1998). Note that the dimensionless cylindrical radius r is scaled to the 
gravitational, rather than Schwarzschild radius such that r ~ R/Rg = Rc^/GM. The four physical disk quantities 
above are determined by four conservation laws. Namely, conservation of mass, vertical momentum (hydrostatic 
balance), angular momentum, and energy (radiative equilibrium). 

That is, we assume that it is possible for all of the angular momentum transport and accretion to take place in the 
body of the disk even if the overwhelming majority of the dissipation resulting from the action of accretion takes place 
in the corona. Thus, in order to constrain these quantities in the corona, we arbitrarily assume a value for the coronal 
optical depth Tc. In the radiation pressure dominated limit, hydrostatic balance may be written as 

orFd^GM^C 

mpc R^ RB ^ ' 

where the left hand side is the upward radiative acceleration associated with a flux F^, and the right hand side is the 
force of gravity in the thin disk limit. Defining He as the height above the midplane where Fd = Q we find 



He 3 . ,V 
— = — nir — . 
R 2 C 

Following SZ9 4, we estimate a coronal density from the scale height and optical depth Tc = Hf.pc'^T /^p- 
Inserting eq. IA6I for He in the radiation pressure dominated case yields 



(A6) 



Pc= —m Tc — , A7) 

For Tc < 1, the radiation pressure can be approximated as 

Prad c ~ A— = -f-mr (A8) 

where A is a factor of order unity which is determined by the angular dependence of the radiation field. 
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B: ESTIMATE OF DISK ALBEDO 

The disk albedo ad reaches large values when e, the ratio of absorption to total opacity, is small. In order to keep 
the photo-electric absorption rate at a low value, the disk must be close to a fully ionized state, which is roughly 
determined by the density of the gas, number of ionizing photons, and recombination rate. A comparison with the 
results of other works is facilitated by the introduction of the "ionization parameter" ^7, given by 

^i^ — ^ . ^ 2x106 — — — [—L] M- - ergcms-i. (Bl) 



jiH Ppd \ P J vo.iy \nbj \ 0.9 J vo-9/ V12 

Here, /3 parameterizes the density of the X-ray absorbing or reflecting region near the surface of the disk and is always 
chosen such that (3 < 1. We have eliminated the relevant relativistic and no-torque correction factors, which take on a 
value ~ 0.003 at r ~ 12. Detailed spectral and ionization calculations show that ad ^ 1 when ^/ ~ 10^ (see Ballantyne 
2002 for a concise review). 

In the following, we motivate why ~ 1 for ^/ ~ 10^ erg cm s^l. For the sake of simplicity, we group all metals into 
a single hydrogenic species with abundance Az- This approximation is only valid when the absorber/reflector is close 
to being fully ionized and when the incident X-ray spectra is roughly described as having equal power across photon 
energy since different ions have different ionization potentials. With this, the expression that governs the balance 
between the number of hydrogenic and completely stripped ions is given by 

-^a {v) dv = neUi+iaBXT) (B2) 

where n^, rii+i, rig, F^, ^nd aii{T) is the number density of bound hydrogenic ions, number density of com- 

pletely stripped ions, number density of free electrons, incident flux, absorption cross section, and recombination rate, 
respectively. By assuming a flat F = 2 spectral index for the ionizing continuum and that (J {v) = ctq {v/vq)^^ an 
expression for the hydrogenic fraction Xi reads 

rii TiH Id h vq (1 ) rie 



Here Fj = Fy and if we further approximate that F ^ Fj then by use of cq. (|B4|I 

for typical values of the micro-physical parameters i.e., ctq ~ 10~^*cm^, hi/Q ^ lOkeV, and afi{T) ~ 10"^" s~^. Now, 
we are in a position to calculate e the ratio of absorption to scattering opacity in terms of the fiducial metal abundance 
Az 

Kabs Xi^ZCTo r:..^n-if \ f 2 X W 



Ksc CTT VSXIO-VV i J 

Note that the value of Az corresponds to Fe at solar abundances. Due to the large disk temperatures in XRBs, the 
relatively abundant CNO metals are fully ionized. Finally, in order to calculate the a^, the disk albedo, we make use 
of the two-stream approximation - an accurate approximation in the elastic limit for photon energies under 10 keV - 
for the disk surface layers (see e.g. HM93). 

1 - (^''^ 1 - 0-024 

ad — 7777 — — 0.953 (B6) 

H-ei/2 1 + 0.024 ^ ' 

for our choice of parameters at r 12. 
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